Cd- Vacancy and Cd-Interstitial Complexes in Si and Ge 
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The electrical field gradient (EFG), measured e.g. in perturbed angular correlation (PAC) experi- 
ments, gives particularly useful information about the interaction of probe atoms like '^^'^In/ '^^'^Cd 
with other defects. The interpretation of the EFG is, however, a difficult task. This paper aims at 
understanding the interaction of Cd impurities with vacancies and interstitials in Si and Ge, which 
represents a controversial issue. We apply two complementary ab initio methods in the framework 
of density functional theory (DFT), (i) the all electron Korringa-Kohn-Rostoker (KKR) Greenfunc- 
tion method and (ii) the Pseudopotential-Plane-Wave (PPW) method, to search for the correct local 
geometry. Surprisingly we find that both in Si and Ge the substitutional Cd-vacancy complex is 
unstable and relaxes to a split-vacancy complex with the Cd on the bond-center site. This complex 
has a very small EFG, allowing a unique assignment of the small measured EFGs of 54MHz in Ge 
and 28MHz in Si. Also, for the Cd-selfinterstitial complex we obtain a highly symmetrical split 
configuration with large EFGs, being in reasonable agreement with experiments. 



I. INTRODUCTION 

The detailed understanding of the electronic structure 
of intrinsic defects in semiconductors is in many cases 
still an unresolved problem. Many different experimental 
methods aiming at microscopic information about intrin- 
sic defect structures exist. The "state of the art" experi- 
ments for these kind of problems are the measurements of 
hyperfine interactions. One example is the perturbed an- 
gular correlation measurement (PAC) of the electric field 
gradient (EFG)^'^i^^i'^. The EFG describes the interac- 
tion of the electrons with the electric quadrupol moment 
Qij of the nuclei and is given by the second derivative 
Vij of the Coulomb potentials at the nuclear site. 

Since the PAC technique itself is based on a radioac- 
tive decay, it is essential for Si and Ge host crystals to 
implant a proper probe atom. The transformation of 
implanted ^^-'^In isotopes to ^^""^Cd by electron capture 
leaves the ^"'^^Cd probe in an excited state and the sub- 
sequent decay to its ground state makes the determina- 
tion of the EFG possible. That means the EFG is just 
measured at the Cd site and reflects the local electronic 
structure which again depends on the atomic environ- 
ment. 

Since the EFG vanishes for cubic or tetrahedral point 
symmetry, probe atoms like Cd give no signal on substi- 
tutional sites. An EFG only arises due to the symme- 
try lowering induced by other defects. In addition the 
direction of the main axis of the EFG gives very clear 
information about the symmetry of the corresponding 
complexes formed. Therefore PAC measurements yield 
particularly useful information about the interaction of 
Cd probe atoms with other defects. The only disadvan- 
tage of the method is that the resulting EFGs are diffi- 
cult to understand and that reliable information about 
the EFG and the electronic and geometrical structure of 
complexes can only be obtained by ab-initio calculations. 

The paper aims at understanding PAC experiments 
with Cd probes in Si2*iaiii and Go^^i^^i". Of particu- 
lar interest is the interaction with intrinsic defects, i.e. 



vacancies and selfinterstitials. In previous studies two 
complexes in Ge have been identified, both with main 
axis along (111) and asymmetry parameter 77 = but 
with very different quadrupolar coupling constants v, i.e. 
1^1 = 54MIIz and 1^2 = 415MIIz. Both values, which dif- 
fer by an order of magnitude, were measured^'^ by two 
different groups and interpreted as Cd complexes with 
intrinsic defects. However the two groups disagree about 
the assignment of the EFG. For Feuser et al. 1^2 is as- 
signed to a substitutional Cd-vacancy complex, while i^i 
should refer to an interstitial Cd atom in the vicinity of 
a vacancy. In contrast to this Sielemann et al.^^ assign 
vi to the substitutional Cd atom with a vacancy on the 
nearest-neighbor site and 1/2 to a substitutional Cd with 
a selflnterstitial on a (111) oriented interstitial position. 
Compared to Ge, the interpretation of the measured fre- 
quencies is much more uncertain in Siii. 

In this paper we perform ab-initio calculations to solve 
this problem. Basically we can confirm the assignment 



of Sielemann et. 



I that the small frequency vi be- 



longs to the Cd-vacancy complex and the high frequency 
to the interstitial one. However in both cases the config- 
urations are very unusual: The Cd atom sits on a bond- 
center position and the single vacancy as well as the inter- 
stitial are split on two neighboring sites. The existence of 
such an impurity-split vacancy configuration has already 
been suggested by Watkinsi^ for a Sn-vacancy pair in 
Si. In a forth-coming paper we will show that this pair 
configuration is also the most stable one for the vacancy 
complexes with other impurities like In, Sn and Sb in Si 
and Ge. Therefore we expect that our results are typical 
for oversized impurities in group IV semiconductors. 



II. THEORETICAL METHODS 

All calculations are based on density functional the- 
ory in the local density approximationi^. Two differ- 
ent methods have been used to solve the Kohn-Sham 
equations. The first one is the pseudopotential plane 
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Table I: Parameters used to construct the pseudopotential projectors. For the KB-potentials only the cut-off radii are listed 
for each 1-channel, for the PAW-potentials the cut-off radii / reference energies (in Ry) are listed. The first line for Cd refers to 
the eigenvalues of the bound states for the respective 1-channel, the second line to the specific energies for the unbound states 
used to construct the second projector for each 1-channel. All radii are in as (0.529177 A). 



wave (PPW) method, which has mostly been used to 
investigate the configuration and stability of the Cd- 
split vacancy and the Cd-split interstitial configurations. 
These configurations, as well as the relaxations of the 
neighboring atoms, have been recalculated by the KKR- 
Greenfunction methodiSiiLiS, which as an all electron 
method allows also to calculate electric field gradients 
and hyperfine fieldsS^iSiSS. 

The PPW method approximates the inhomogeneous 
system containing defect complexes by periodically ar- 
ranged large supercells and uses plane waves to ex- 
pand the electronic wave functions. This has the 
advantage that band-structure methods can be used 
to determine the electronic structure, and total en- 
ergies and forces on the atoms can be calculated 
without difficulty for arbitrary arrangements of the 
atoms in the supercell. Our EStCoMPP-program^^ 
can handle norm-conserving Kleinman-Bylander (KB)- 
pseudopotentials??' and projector-augmented (PAW)- 
multi-projector-pseudopotentials^'^ used to describe the 
electron-ion-interaction. For Si and Ge KB-non-local 
pseudopotentials are used for the s,p- valence electrons 
and the d{l = 2)-component is used as a local poten- 
tial. For Cd the 5s, p- valence electrons as well as the the 
4d-electrons are treated explicitly with PAW-potentials 
and the f{l = 3) -component is used as a local potential. 
Partial-core-correction was used for Gd. The atomic con- 
figurations, cut-off radii and reference energies used to 
generate the projectors are listed in Table Q] For the ac- 
tual calculations we used a (lll)-oriented supercell with 
the basis vectors hi — 3a(0, —1, 1); 62 — 3a(— 1, 1, 0); 63 = 
2a(l, 1, 1) containing 108 atoms, a is the theoretical lat- 
tice constant ( 10.44 as for Ge and 10.25 a^ for Si; Bohr 
radius as — 0.529177 A). The Cd-complex was placed in 
the middle of the cell. We used Cs^ symmetry explicitly 
for all configurations. A plane- wave basis set equivalent 
to 6 X 6 X 6 Monkhost-Pack^'' fc-points was used which 
yields three inequivalent fc-points in the irreducible part 
of the Brillouin-zonc for our supercell. We used a plane- 
wave cut-off of 20.25 Ry, which is required for the con- 
vergence of the Gd PAW-potential. The atoms belonging 
to the Gd-complexes and their nearest neighbors were 
relaxed until the forces on all atoms were less than 0.1 
mRy/as- We checked that forces on further neighbors 
(which were not moved) were of the order of 1 mRy/as 



and we estimate that the positions of all relaxed atoms 
were determined with an accuracy better than 0.01 as- 

In the KKR Greenfunction method the calculation is 
devided in two steps. First the Greenfunction of the host 
is determined. In a second step the host Greenfunction is 
used to determine the Greenfunction of the crystal with a 
single impurity by a Dyson equation. For details we refer 
to referenceSSr. All calculations are performed with an an- 
gular momentum cut-off of Imax = 4. For the host crystal 
we use the LDA lattice constantsS& of Si (10.210^) and 
Ge (10.53aB). A k-mesh of 30 x 30 x 30 fc-points in the 
full Brillouin zone is used. We use the screened KKR- 
formalism2L2&, with the TB-structure constants deter- 
mined from a cluster of 65 repulsive potentials of 4 Ry 
heights. The diamond structure is described by a unit cell 
with 4 basis sites, two for host atoms and two for vacant 
sites. For the Greenfunction of the defective systems, we 
allow 77 potentials of the defect and the surrounding host 
atoms to be perturbed, which are then calculated selfcon- 
sistently with proper embedding into the host crystal. All 
calculations include the fully anistropic potentials in each 
cell and thus allow the reliable calculation of forces, lat- 
tice relaxations and electric field gradients. The single 
particle electron density n{r) can be calculated via the 
imaginary part of the Greenfunction: 
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The Goulomb potential can be expressed in terms of n{f) 
and therefore the force F" on atom n can be derived with 
the use of the "ionic" Hellmann-Feynman theoremS^: 
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where the first term denotes the Goulomb force of the 
valence electrons and of the other nuclei on the nuclear 
charge Z", while the second term gives the force on the 
core electrons (G) of atom n, evaluated by the derivative 
of the Kohn-Sham potential. It can be shown, that in an 
angular momentum expansion the force is determined via 
the / — 1 component of the charge density. The tensor of 
the EFG is given by the second derivative of the Goulomb 



potential: 
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a,P^x,y,z; Tr V^/j = 0, 



where (x, y, z) denotes the coordinate system. In practi- 
cal applications, since we expand all quantities in spher- 
ical harmonics, the electric field gradient is calculated as 
followsi2*Sa: 
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Figure 1: The Cd-split vacancy configuration 
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where l2m(?0 are spherical harmonics, V2,mir) are the 
I = 2 components of the Coulomb potential and Rmt is 
the Muffin Tin Radius of the probe atom. Furthermore 
n2m{'f) are the I = 2 components of the radial charge 
density. We denote that just the I — 2 components are 
needed for the EFGs, while the I = 1 components enter 
into the forces. Due to this, the EFG vanishes for sites 
with cubic or tetrahedral symmetry. For the cell division 
of the crystal we used a generalized Voronoi construction. 
The impurity cell at the bond center was constructed 
such that it is slightly larger than the cells of the host 
atoms, thus avoiding to decrease the Muffin Tin Radius 
of real atoms by more than 16%. For the energy con- 
tour in equation (Q all d-states are included. For the 
treatment of occupied states in the band gap, we have 
used a separate energy contour for gap energies. Since 
we used group theory we are able to occupy these gap- 
states according to their symmetry, in this way avoiding 
to calculate the wave functions of the localized gap states. 



calculations show that this configuration is unstable and 
does not represent a local minimum. Thus a comparison 
of the EFGs with the measured frequency 1/2 = 415 MHz 
in Ge is meaningless, since the configuration is not stable. 
Apparently the instability arises from the large size of the 
Cd atom, which is pushed away from the three Si or Ge 
neighbors into the empty space towards the vacancy and 
settles on the symmetrical bond-center position. 

The divacancy in Si and Ge: Since the electronic struc- 
ture of the Cd-split vacancy complex can most easily be 
understood by first considering the divacancy, we present 
in Fig|21the calculated density of states (DOS) of the di- 
vacancy in Si. Four different peaks can be seen: single 
degenerate Aig and A2« resonances in the valence band, 
and double degenerate Eg and Eu localized states in the 
gap. The splitting into g and u state arises from the 
hybridization of the A- and E-dangling bond states of 
the two single vacancies into bonding and anti-bonding 
combinations. In a molecular orbital (MO) scheme the 
energy levels are schematically indicated in Fig|31on the 
right. In the calculation the Fermi level is situated at the 
maximum of the valence band E^ai ■ Therefore only 4 of 
the 6 dangling bond electrons can be accommodated, so 
that in the calculation the charge state of the divacancy 
is2-h. 



III. CD-SPLIT VACANCY COMPLEX 

The substitutional Cd-vacancy complex: The starting 
point of our study is the substitutional Cd-vacancy com- 
plex. To our surprise this complex with the Cd on the 
substitutional site and the vacancy on the nearest neigh- 
bor site is highly unstable. When fixing the Cd at the 
substitutional position and allowing the six host atoms 
to relax, the energy of the configuration is for both Si 
and Ge hosts more than 1 eV higher than the energy 
of the stable split configuration shown in Fig^ When 
the Cd atom is allowed to relax, the energy continuously 
lowers (without energy barrier) until the Cd reaches the 
stable bond-center position shown in Fig^ The fixed 
substitutional configuration is highly anisotropic and the 
calculations give rather high EFGs for the Cd atom, with 
values of about -420 to -430 MHz in Ge and -390 to -420 
MHz in Si, where the spread is caused by the different 
charge states. However both the PPW and the KKR 
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Figure 2; Local density of states (LDOS) for the divacancy'^"'" 
in Si. Shown are the local s, p and d contributions on a 
vacancy site, separately for the Aig and A2u subspaces(left) 
and the Eu and Eg subspaces(right). As can be seen, for the 
A subspace we find mainly s-contributions and for the E sub- 
space mainly p-contributions. The other 1-channels contribute 
less than 10% over the entire energy range. 
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Cd-atom complex divacancy 

Figure 3: Schematic tight-binding scheme of energy levels to 
explain the electronic structure of the Cd-split vacancy com- 
plex in Si and Ge. Included are the 5s, 5p levels of the Cd 
atom and the dangling bond states of the divancy in its 2+ 
state. From left to right: the single Cd-atom, the Cd-split 
vacancy complex and the divacancy^"*" 




Figure 4: LDOS projected on the irreducible subspaces of the 
Dsd group at the Cd site (be) for the Cd-split vacancy in Si 
(above) and Ge (below). The nearest neighbors are fixed on 
ideal lattice positions. 



Cd-split vacancy complex: The local DOS at the Cd 
bond-center site of the Cd-split vacancy complex is shown 
in Fig^ both for Si and Ge. Since this complex preserves 
the point symmetry D^d of the divacancy and since the 
Cd 5s and 5p states are well above the Fermi level, the 
decomposition and level structure of the divacancy DOS 
are basically preserved. As indicated by the MO level 
scheme in the center of FigEl , by hybridization between 
the states of the Cd and of the divacancy the divacancy 
levels are pushed to lower energies and the Cd states to 
higher energies. The selfconsistent calculations yield the 
Eu state as a resonance at the top of the valence band, 
in agreement with the level scheme assumed for the split 
complex in FigOl In this configuration 8 electrons are ac- 
commodated in the single degenerate Aig and A2u states 
and the double degenerate E^ state. These 8 electrons 
just correspond to the 6 electrons of the dangling bonds 
of the divacancy and the two electrons of Cd, so that the 
calculations yields the neutral charge state, as long as the 
Eg level in the gap remains unoccupied. The d-states of 
Cd are located around —8.5 eV and, while included in 
the calculations, are basically localized and of minor im- 
portance for the electronic structure and the EFG. 

As Fig^shows, the Cd atom has in its nearest coordi- 
nation shell six host atoms at a distance at 1.26 times the 
nearest neighbor distance (if lattice relaxations are ne- 
glected; see below) . Therefore the hybridization between 
the Cd 5sp states and the divacancy states is rather weak, 
so that the occupied Aig, A2u and Eu states contain lo- 



cally only relatively small admixtures of these Cd states, 
leading to a charge of only 1.0 electrons (in the case of 
Ge; for Si: 0.9 electrons) in the Cd sphere, which has 
about the same size as the Si or Ge spheres. Therefore 
we can understand the electronic structure in a simple 
picture in real space, by considering the Cd as a "naked" 
Cd'^^ ion embedded at the bond-center position into the 
divacancy. The strongly attractive Cd^+ potential low- 
ers all divacancy states, but preferentially, and in first 
order perturbation theory only, the Aig level which has 
the full D^d symmetry of the defect. Due to this, the Aig 
state is stronger localized at the Cd site than the other 
states and shows up more prominently in the LDOS of 
Cd. The LDOS for the Si and Ge hosts are qualitatively 
and quantitatively very similar. 

The different charge states [Cd Y]^' , [Cd V]^-, ... of 
the Cd-split vacancy complex are obtained by occupying 
the double degenerate Eg state in the gap. The [Cd V]^^ 
state is magnetic and exhibits a small moment about 0.02 
fiB at the Cd site, both in Si and Ge. Since the total 
moment is 1 /is, this shows the large spatial extend of 
the Eg state. 

Lattice Relaxations: FigEl shows the calculated lattice 
relaxations of the six nearest neighbor Si-atoms around 
the Cd atoms, while the table gives the calculated values 
for the case of Ge. All relaxations are given in percent- 
ages of the nearest neighbor distance of the host crystal. 
First, we see that in the case of the neutral complexes 




5.2% 



6.2% 



Defect 


KKR [% NN] 


PPW [% NN] 


\V\Cd\Vf in Ge 


6.6 


6.0 


[V|Cd|V]'- in Ge 


7.5 




\V\Cd\V\'- in Ge 


8.5 





Figure 5: Relaxations of the 6 nearest neighbors of the Cd- 
atom, Cd-split vacancy complex in the figure for the Si host, 
in the table for the Ge host. Besides the KKR values for the 
neutral configuration also the PPW values are given. 



the KKR-values agree well with the PPW ones. The 
small differences arise from the different ways used to de- 
scribe the electron-ion interaction (pseudopotentials for 
the PPW and all-electron in KKR) and the shghtly differ- 
ent lattice constants. They are of the same order of mag- 
nitude as typical differences in other instances, e.g. for 
the relaxations of the near surface atoms of Cu surfaces. 
Secondly, we notice that the relaxations point toward the 
Cd atom and increase by about 1 % with increasing total 
charge of the complex. Both facts are consequences of 
the electrostatic interaction, i.e. the Cd^+ ion attracts 
the dangling bond electrons of the Si or Ge neighbors 
and the attraction increases with increasing charge in 
the extended Eg states. Third, the relaxations in Ge are 
larger than the relaxations in Si, which is a consequence 
of the 3 % larger lattice constant of Ge as compared to Si 
(theoretical values), so that the local distances are very 
similar in both hosts. Despite the sizeable relaxations 
the final distances are still very large, the strength of 
the hybridization is still weak and the LDOSs are hardly 
changed by relaxations. 

Charge anisotropy and EFGs: As stressed in the intro- 
duction the EFG is determined by the anisotropy of the 
charge density at the Cd site, more correctly by the £ — 2 
components of the (^r7T,)-projected charge density n2m{r)- 
Such £ = 2 components arise from the anisotropy of the 
p- and the d-charge densities, since mixed s — d ot p — f 



terms are very small. In a good approximation the EFG 
is therefore determined by the anisotropy of the local p 
and d chargaiSiSS 



T4, = aAN^ + bANd 



(5) 



where ANp and AA^^ are the anisotropic p and d charges, 
which for the considered symmetry and z-axis in (111) 
direction are given by 

AN,^ J.^'^dE'll {n,^{E') + n,^{E'))-n,AE')] 
AN, = /„^- dE' K^,_^, iE') + na^^ (E') (6) 



The positive constants a and b in equation JHJ are matrix 
elements of -^ with normalized p and d wave functions 
at Ep- 

We can now analyze the anisotropic charges in terms 
of the contributions from the different irreducible contri- 



butions Alg, A-), 



E,, 



and Eg to the local Cd DOS as 



^xy-j ^: 



yi and dxy are very small, so that 
and dyz are 



given in FigQl The Aig subspace only contains the local 
s and dz'i contributions, the subspace A2u only the pz, 
the i^ti-subspace the Px and py and the Eg state the local 
j,2, dxz and dyz contributions. 

In the neutral state [CdV]*^ the Eg state is not occu- 
pied and we find negative electric field gradients. The 
contribution ANp is negative, since the £"„ state cannot 
overcome the negative Up^ contribution of A2u and AiV^ 
is negative due to the n, ^ contribution from Aig. When 
one then populates the Eg-siaXe, we find that the con- 
tributions from d^-i 

additional negative contributions from d. 
expected for the charged configurations 

The calculated values for the relaxed configurations are 
given in tabled First, we would like to point out that 
the displacements of the NN-atoms play only a minor 
role, which is not surprising in view of the small changes 
found in the LDOS. Secondly, we observe that for the 
charged defects the EFGs in Si and Ge increase for each 
electron added in the Eg state by about 20 MHz, more 
correctly about 21 in Si and 23 in Ge. This is a con- 
sequence of the fact that the Eg state is very extended: 
For each electron added to the complex the local charge 
of Cd increases only by 0.05 electrons, which represents 
a weak perturbation. Finally we find that the EFG of 



host crystal 


Si 
EFG [MHz] 


Ge 

EFG [MHz] 


Si 
HF [kG] 


Ge 

HF [kG] 


[V\Cd\Vf 


-6.97 


-32.69 






[V\Cd\V]'' 


-27.99 


-55.69 


21.8 


8.5 


[V\Cd\Vf- 


-49.47 


-79.64 






EXP 


± 28.00 


± 54.00 







Table II: Calculated EFG on fully relaxed positions for the 
Cd-split vacancy in Si and Ge for different charge states 
(0, 1 — ,2— ). We used Q=0.83barn for Cd. Experimental val- 
ues (EXP) are taken from reference^*-. 



the singly charged configuration [CdT^]^" agrees both in 
Si and Ge very well with the measured values. Note 
that experimentally the sign of the EFG cannot be de- 
termined. We consider the close agreement as a confir- 
mation that the Cd-split vacancy complex is the most 
stable Cd-vacancy complex and that the measurements 
refer to the negatively charged [CdV]^~ complex. This 
is in agreement with Haesslein et ali, who assigned the 
frequency vi = ±54 MHz in Ge to a Cd-vacancy complex 
in the charge state 1 — . As mentioned earlier already, in 
this configuration the complex is magnetic and for the 
Cd-atom our calculations predict a small hypcrfine field 
(HF) of 21.8 kG in Si and 8.5 kG in Ge. Up to now these 
hyperfine fields have not been measured. In summary we 
find that the small EFGs measured are a typical feature 
of the Cd-split vacancy configuration, in which the Cd 
atom is far away from the nearest neighbors so that the 
hybridization is very small, much smaller than for the 
strongly anisotropic (but unstable) substitutional config- 
uration, yielding EFGs of around 400 MHz. Partly the 
geometrical structure also helps, since the Cd atoms and 
the six neighbors form a distorted octahedron. For an 
ideal octahedron the EFG would vanish due to symme- 
try. 



IV. CD-INTERSTITIAL COMPLEX 

The Cd-selfinterstitial complex is a possible candidate 
for the large EFG of 1^2 = 415 MHz measured by two 
groupsiiSi in Ge. In our electronic structure calculations 
we have first studied the interstitial-Cd complex shown 
in FigEl (upper left): a self-interstitial on a tetrahedral 
site adjacent to a substitutional Cd impurity. To our 
surprise we find again, that this configuration is very un- 
stable, exhibiting a large force of about 100 mRyd/as 
on the Cd atom. By relaxing this configuartion the three 
adjacent atoms in (111) direction, i.e. the selfinterstitial, 
the Cd atom and the nn host atom in (111) direction 
move more or less uniformly until the Cd atom settles in 
the bond-center position and the two host atoms about 
halfway between the interstitial and the substitutional 
positions. Therefore we refer to this as the Cd-split in- 
terstitial configuration, which is highly symmetric {D^d, 
symmetry). It is shown in FigEl (upper right), together 
with the exact positions of the neighboring atoms (FigEl 
(bottom)). 

The electronic structure of this complex can again be 
understood by considering the Cd atom as a Cd^+ ion. 
FigOshows the LDOS in the cell of the Cd-atom for this 
configuration, in FigC^ without the presence of the Cd 
atom and in Fig[7|D in the presence of the Cd. Only the 
LDOS in the Aig and A2U sub-spaces are shown, since 
the local Eu and Eg contributions are small and not of 
relevance for the EFG. These states are mostly localized 
at the two interstitials and the 6 neighboring host atoms. 

For the two shifted host atoms without Cd, we see in 



FigCK two sharp peaks oi Aig and A2U symmetry at Ep. 
These are symmetric and antisymmetric combinations of 
Pz-dangling bond states of the two split atoms, point- 
ing into the interstitial region. A contour plot in the x-z 
plane of the DOS of these states is shown in Fig|SJi. Since 
the two dangling bond states do practically not overlap, 
the Aig and A2U states are nearly degenerate. By in- 
serting the Cd^+ ion, the bonding Aig state is shifted 
to lower energies and becomes more localized at the Cd 
site, while the A^ states is only slightly shifted to lower 
energies (Figda). Both the Aig and A2U states are fully 
occupied yielding a neutral configuration. The contour 
plot of the A-state close to Ep, being dominated by the 
v42M-state, shows Cd-pz orbitals forming a bridge between 
the two dangling bond states, with a nodal plane perpen- 
dicular to the z-axis (FiglSJ)). 

The calculated EFG of the Cd-split interstitial com- 
plex is quite large, i.e. -395 MHz in Ge and -345 MHz in 
Si. The strong anisotropy arises from the two split inter- 
stitial atoms. The major contribution to the EFG arises 
from the Pz orbital in the ^2u-state, shown in the contour 
plot of FiglHb; in addition small contributions come from 
Px and Py orbitals belonging to i?„ states. The Aig state 
gives nearly no contribution, since the dominant s-state 
is isotropic and the ^^2 -admixture is small. The calcu- 
lated EFG of -395 MHz for the Cd-spit interstitial in 
Ge is in good agreement with the observed frequencji^ 
1/2 = ±415 MHz, which has been assigned by Haesslein 
et alii to a Cd-interstitial, but by Feuser et alr^ to a Cd- 
vacancy complex. The calculated value of —345 MHz for 




56%NN 
(57%NN) 

3.4%NN 
(4 .7%NN) 

Figure 6: The Cd-selfinterstitial complex with Cd on the sub- 
stitutional site (above left) leads after relaxation to the sym- 
metrical complex with Cd on the bond center and the two 
host atoms shifted half way between the substitutional and 
interstitial positions (above right). The lower figure gives the 
calculated displacements of the neighbors from the ideal po- 
sitions in percent of the NN distance for the Ge host. In 
brackets the results of the PPW calculation are given. 
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Figure 7: LDOS in the bond-center cell for the cases that it 
is unoccupied by the Cd impurity (upper figure, a) and that 
it is occupied by Cd (lower figure, b). The position of all 
other atoms are the same in both figures as given for the split 
configuration in Fig. 6 (upper right). 



the interstitial complex in Si is more difficult to compare 
with the experimental data, which have been reviewed by 
Wichert et ali^i. From our calculations we can exclude 
the axial symmetric and < 111 > oriented complex De3 
with ±448 MHz; moreover we can also exclude the de- 
fects De2 since it is not axial symmetric, and Del since it 
refers to the Cd-split vacancy complex. A good match is 
only obtained with the unidentified defect X2 with ±334 
MHz, being axial symmetric and < 111 >-oriented. 



V. SUMMARY 

In this paper we have investigated the electronic and 
geometrical structure of complexes of the radioactive 
probe atom Cd with intrinsic defects in Si and Ge. In 
particular we have calculated the electric field gradient 
(EFG) of the Cd probe atom in these complexes, with the 
aim of assigning measured EFG frequencies to calculated 
defect complexes. The calculations are based on density 
functional theory in the local density approximation and 
apply the full-potential KKR Green's function method 
as well as the pseudo potential plane wave method. Both 
methods give nearly identical results for the atomic con- 
figurations. 

The most simple complex, a substitutional Cd-atom 
adjacent to a vacancy on a nearest neighbor site, is found 



to be unstable and to relax into a Cd-split vacancy config- 
uration: the Cd-atom sits on the bond centered site, with 
the vacancy split into two half- vacancies on the neighbor- 
ing substitutional sites. As seen from the Cd atom, this 
complex is very isotropic and results in a rather small 
EFG. The values for the single negatively charged com- 
plex, i.e. —55.7 MHz in Ge and —28 MHz in Si, are in 
excellent agreement with measured EFGs (±54 MHz in 
Ge and ±28 MHz in Si) and also the charge state agrees 
with conclusions of Haesslein et ali. This close agree- 
ment confirms the existence of the calculated Cd-split 
vacancy configuration with very small EFG. By investi- 
gating the most simple complex of a Cd-atom with an 
interstitial, i.e. a substitutional Cd atom with a selfin- 
terstitial on the adjacent tetrahedral site, we find that 
this complex is also unstable and decays spontaneously 
into a Cd atom on the bond-center position and the two 
host atoms shifted from the substitutional site halfway 
to the tetrahedral position, thus forming a Cd-split in- 
terstitial complex aligned in < 111 >-direction. For Ge 
the calculated EFGs of —415 MHz agrees well with the 
experimental value of ±395 MHz assigned by Haesslein 
et al. to a Cd-interstitial complex. The experimental 
situation in Si is not so clear. 

These unusual split configurations are typical for in- 
trinsic complexes with oversized impurities, as we will 
demonstrate in a future paper for a series of vacancy 




Figure 8: (color-online) A contour plot in the xz-plane of the 
A2u state close to the Fermi energy (see FigUJ for the Cd split 
interstitial complex: Fig. 8a (upper panel) in the absence of 
the Cd impurity. Fig. 8b (lower panel) in the presence of the 
Cd on the bond-center position. 



complexes with oversized and "normal" impurities. 



for suggesting the Cd-split vacancy configuration to us. 
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